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Abstract. 

The beginning of the cosmological phase bearing the direct kinematic imprints 
of supernovae dimming may significantly vary within different models of late-time 
cosmology, even if such models are able to fit present SNe data at a comparable level 
of statistical accuracy. This effect - useful in principle to discriminate among different 
physical interpretations of the luminosity-redshift relation - is illustrated here with a 
pedagogical example based on the LTB geometry. 

It is by now widely known that the observed luminosity-redshift distribution of 
type la supernovae [1] can be fitted even without dark energy, provided one introduces a 
sufficiently inhomogenous space-time geometry. A typical, very simple example of such 
a possibility is provided by matter-dominated cosmological models of the Lemaitre- 
Tolman-Bondi (LTB) type (see e.g. [2] for an incomplete list of papers on this 
subject), provided the observer is located near enough to the symmetry centre of the 
inhomogeneous - but isotropic - matter distribution [3, 4]. 

Such an example may be regarded as unnatural because of the amount of fine tuning 
required to localize the observer position [4] , and also appears theoretically disfavoured 
by the possible presence of weak geometric singularities [5]. Nevertheless, the possible 
role of inhomogeneities in determining (or at least substantially contributing to) the 
large-scale dynamics should be - and, indeed, currently is [6] - seriously scrutinized and 
discussed, even in the presence of a dominant dark-energy cosmic component [7]. 

The general question that arises in such a context is how to distinguish different 
successful analyses of the SNe data based on different physical models and, in 
particular, on different (homogeneous versus inhomogeneous) large scale geometries. 
For inhomogeneous models of the LTB type various answers to this question are known, 
concerning the size of the voids described by the LTB geometry [8], the local expansion 
rate inside the voids [9], and the associated effect of redshift drift [10]. Other possibilities 
to test LTB models are provided by studies of scalar perturbations [11], of small-scale 
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CMB effects [12], of the cosmic age parameter [13], and of BAO (baryon acoustic 
oscillations) data [14]. 

The main purpose of this note is to point out another possible difference between 
inhomogeneous and more conventional interpretations of the SNe data, not yet discussed 
in the literature; such a difference is based on the value of the redshift parameter z acc 
(to be defined below, see after Eq. (14)), marking the beginning of the regime directly 
characterized by the kinematic imprints of SNe dimming. The value of such parameter 
can be largely different even within models able to fit the presently observed luminosity- 
redshift distributions at a comparable level of accuracy (see e.g. [15] for earlier studies 
on the beginning of the accelerated regime in the context of a homogeneous geometry). 

This suggests two possible experimental ways of discriminating among models of 
the luminosity-redshift relation. First, direct observations able to extend our present 
knowledge of the Hubble diagram up to values of z higher than those allowed by present 
SNe data: for instance, gamma-ray burst (as discussed in [16]), or even gravitational 
waves observations, through an analysis of the luminosity distance of the so-called 
"standard sirens" [17]. Second, indirect observations which are sensitive to the time- 
dependence of the so-called "transfer function" [18], which controls the evolution of the 
primordial perturbation spectrum inside the horizon down to the present epoch, and 
which is crucially affected by the kinematics of the cosmological background (see e.g. 
[19])- 

The possible relevance of the parameter z acc will be illustrated in this paper by a 
simple exercise, in which the SNe data of the recent Union2 compilation [20] are fitted 
using a inhomogeneous, matter-dominated LTB model, and such a fit is compared with 
the standard one performed in the context of the flat concordance ACDM model. We 
stress that our aim is not to provide a realistic alternative to the successful concordance 
cosmology, but only to discuss how to distinguish, at least in principle, different fits of 
SNe data based on different geometric schemes. The proposed diagnostic may be added 
to other general methods aiming at discriminating the expansion history of competing 
models, like - in particular - dark-energy based diagnostics for homogeneous models 
[21]; a Friedmann equation diagnostic for homogeneous versus inhomogeneous models 
[22, 23]; and the already mentioned test of redshift drift [10, 23, 24]. 

The cosmological configuration we will consider refers to a late-time (in particular, 
post-reionization) Universe, characterized by a stochastic distribution of many overdense 
and underdense regions, of various possible sizes and shapes, possibly even incoherently 
superimposed among each other*. Let us suppose that in such a context, and up to a 
given scale r v (to be specified below), the effective (averaged) large-scale geometry can 
be locally described by a model of the LTB type. Such a model is characterized in general 
by three arbitrary functions of the radial coordinate (see e.g. [26]). For the illustrative 
purpose of this paper, however, it will be enough to consider a simple example where 
the contribution of the spatial curvature is negligible and the gravitational sources are 

* Such a configuration is in principle different from that of a typical "Swiss cheese" scenario, where 
the void regions are more or less regularly distributed and well disconnected (see e.g. [25]). 



Discriminating different models of luminosity-redshift distribution 



3 



dominated by an isotropic cold dark matter (CDM) distribution (but the model could 
be easily generalized by the addition of an arbitrary cosmological constant). 

We will assume that the large-scale geometry around a given observer is described 
- in polar coordinates and in the synchronous gauge - by the following metric, 

ds 2 = dt 2 - A' 2 (r, t)dr 2 - A 2 (r, t) (d6 2 + sin 2 Od^ 2 ) , (1) 

where a prime denotes partial derivatives with respect to r and a dot with respect to 
t. In the limit A(r,t) = ra(t) one recovers the well known, spatially flat, Friedman- 
Lemaitre-Robertson- Walker (FLRW) metric. In general, the unknown function A(r,t) 
is to be determined by the Einstein equations, which in our case reduce to 

H 2 + 2HF = 8nGp, 2H + 3H 2 = 0, (2) 

where H(r,t) = A/ A and F(r,t) = A' /A'. The density profile of the CDM distribution 
around a central observer, p = p(r, t) satisfies the covariant conservation equation: 

p+(2H + F)p = Q, (3) 

while all the other Einstein equations are identically satisfied by the metric (1) (see e.g. 
[27]). 

The above cosmological equations can be integrated exactly, and in this paper we 
will adopt the particular exact solution 



A(r, t) — r 



1 + \tH Q {r) 



2/3 



(4) 



normalized in such a way that A = r at t = 0. The arbitrary function H (r) depends 
only on the radial coordinate, and the usual matter-dominated FLRW solution is exactly 
recovered in the limit H = const. We will use, in particular, the parametrization 

H Q {r) = H + AHe~ r ^ v , (5) 

already suggested in [27] for a similar LTB scenario (a brief discussion of other possible 
choices for the phenomenological profile H (r) will be given in the final part of this 
paper). For the chosen profile the combination of parameters H + AH = H o (0) 
corresponds to the locally measured value of the Hubble constant, while the distance r v 
represents the typical distance scale above which inhomogeneity effects become rapidly 
negligible. 

To make contact with more general forms of the LTB metric appearing in the 
literature, and expressed in terms of three functions M(r), ts(r), E(r), (we are following 
the notations of [26]), it may be useful to report here the values of those functions for 
the model we are using. The effective gravitational mass with comoving radius r, for 
our solution, is given by M(r) = (l/2)r 3 H^(r). It can be easily checked that this 
function grow like r 3 for r <^ r v and r ^> r v , while, in the transition regime r ~ r v , 
it is characterized by a fractal index D = 0.4, i.e. M{r) ~ r 3 - D r 2S . The time 
scale t B - i.e., the local "big-bang time" at which A(r,t) = - in our case is given by 
t B (r) = -(2/3)H 1 (r). 
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Finally, it is important to stress that the obtained solution is consistent with our 
assumption of vanishing spatial curvature, i.e. with the choice E(r) = 0. Perturbing 
the solution with the addition of scalar curvature (and assuming that E(r) ~ r 2 as in 
the large-scale FLRW limit), we have checked indeed that the curvature contribution 
to the total energy density may have a variation which is at most of the order of 0.05% 
over length scales of order r v and time scales of order Hq 1 . Hence, if initially small but 
nonzero, it keeps small over the whole spatial and temporal range of interest for this 
paper. 

Let us now compute the luminosity distance d^ of a source emitting light at a 
cosmic time t and a radial distance r from the origin. We will assume, for the moment, 
that the observer is also located at the origin (the consequences of a possible off-center 
position will be discussed later). The angular distance (or area distance) of the source, 
for the metric (1), is then given by dA = A(r,t), and the luminosity-distance, according 
to the so-called "reciprocity law" [28], reduces to d L = (1 + z) 2 A(r,t), where z is the 
redshift parameter evaluated along a null radial geodesic connecting the source to the 
origin. 

Calling the static (time-like) geodesic vector field tangent to the worlines of 
source and observer, and /c M the null vector tangent to the null radial geodesic, we find 
in our metric = dx^/dr = (1,0,0,0) and jfe** = ((A'y 1 , -(A'y 2 , 0, 0). Hence, for 
light emitted at time t, radial position r, and observed at the origin at t — t , 

(^) , t0 A '^y 

where A' = A'(0,t ) = const. 

For the phenomenological applications of this paper we need to express d^ 
completely in terms of the redshift, namely we need to invert Eq. (6) to determine 
r(z) and t(z). We may consider, to this purpose, the differential variation of z with 
respect to the proper time interval dr separating two different instants of light emission, 
at fixed observation coordinates: dz/dr = u^d^z = — (1 + z)A'/A'. It follows that, along 
a null radial geodesic (where dt = —A'dr): 

dt _ dtdr _ A' dr I dt 1 

dz~ drdz~ (l + z)A r dz ~ A' dz ~ (l + z)A r 

For the model of Eq. (4), in particular, we obtain the differential equations 

dt 1 [2 + 3tH (r)\ [2 + 3t H (r) + 2rtH' Q (r)] 

dz ~ " 1 + z QtH 2 (r) + ArH' {r) + 4# (r)[l + rtH' (r)} ' 

dr = _^ [2 + 3tH (r)}^ 

dz 2V3(1 + z ) 3tH$(r) + 2rH' (r) + 2H (r)[l + rtH' {r)} ' 1 ' 

Solving the above equations for t(z), r(z), and inserting the solutions into the explicit 
definition of d^, 



d L (z) = (1 + z) 2 A(r, *) = (! + z? r(z) 



1 + h(z)H (r(z)) 



2/3 



(9) 
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fi(z) = 5 log 



10 



+ 25. (10) 



we are now in the position of comparing the predictions of our model with the 
observational data (as well as with the predictions of the standard ACDM scenario). 

Let us first recall that the Union2 compilation of the Supernova Cosmology Project 
[20] concerns redshift-magnitude measurements of 557 SNe of type la and provides, for 
each supernova, the observed distance modulus (with relative error) fi ohs (zi) ± Afj,(zi), 
i — 1, . . . , 557, for redshift values ranging from z\ = 0.015 to z 557 = 1.4. The distance 
modulus fj,(z) controls the difference between apparent and absolute magnitude, and is 
related to the luminosity distance <1l{z) by: 

d L {z) 
lMpc_ 

Here di is given in units of Mpc, and the constant number 25 is determined by the 
conventional reference scale assumed for the absolute magnitude. 

The luminosity distance of Eq. (9), with Hq{t) given by Eq. (5), is characterized in 
principle by three independent parameters, and can be applied to fit the experimental 
data by allowing free variations of H, AH and ry. We have performed that exercise, 
and found that the resulting best fit provides for H (0) = H + AH a value very close 
to 70 Kms _1 Mpc _1 . We have thus chosen to concentrate the present discussion on a 
simpler, two-parameter fit of the data - which, in any case, is sufficiently accurate for the 
illustrative purpose of this paper - by imposing on our model the "a priori" constraint 
H + AH = 70Kms~ 1 Mpc~ 1 . In this way we can eliminate, for instance, H, and we can 
fit the experimental points fi ohs (zi) ± A/i(zi) by performing a standard x 2 analysis with 



557 



2 



v2 = v- /i obs (^)-M^rv,Aijy 
x UY am*) J ' 1 } 

The theoretical values fx(zi, ry, AH) can be determined, for each value of Zi, by 
numerically integrating the two equations (8), and computing the corresponding d^Zi) 
as a function of the two parameters ry , AH. By minimizing the above \ 2 expression we 
have found the best fit values 

r v = 3000 ±497 Mpc, AH = 26.6 ± l.SKms^Mpc" 1 , (12) 

at a confidence level of 95%, and with a goodness of fit x 2 /d.o.f. = 0.99. The 
minimization has been performed using the MINUIT package from CERNLIB [29]. The 
result of the fit is graphically illustrated by the red curve plotted in the left panel of 
Fig. 1, superimposed to the full set of Union2 data (reported with error bars). 

Consider now, for comparison, a fit of the same data performed in the context of 
a spatially flat FLRW geometry, with perfect fluid sources representing CDM and a 
cosmological constant A. Denoting with Q m and Q A the present critical fraction of dark 
matter and dark energy, we can express the luminosity distance in the usual integral 
form as 

d L (z) = ±±1 f dx [n rn (i + x f + n A ] - 1/2 (is) 

-no Jo 

(see e.g. [30]). Proceeding as in the previous case, we will reduce the number of 
parameters from 3 to 2 by imposing the same phenomenological constraint as before, 
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which in this case amounts to the condition Ho(Q m + S^a) 1 / 2 = 70 Kms _1 Mpc _1 . Using 
Eq. (13) to compute fi(zi, fl m , fl\), and minimizing the corresponding \ 2 expression, we 
obtain the best fit values Q m = 0.27 ± 0.01, ^a = 0.71 ± 0.03, at a confidence level of 
95%, with x 2 /d.o.f. = 0.98. The result of the fit is illustrated by the blue curve on the 
right panel of Fig. 1. 




0.5 1 1.5 0.5 1 1.5 



z z 

Figure 1. The Hubble diagram of the Union2 dataset. The left panel illustrates the 
best-fit result for a two-parameter fit of our example of inhomogeneous geometry, with 
x|, TB /d.o.f. = 0.99. In the right panel we present the corresponding best-fit result for 
a homogeneous ACDM model, with Xacdm/^- ^ ■ = 0.98. 



The luminosity-redshift relations of the two models of Fig. 1 are in good agreement 
with the data, and in both cases the data points are fitted at a comparable level 
of statistical accuracy. However, we can disclose an important physical difference 
between the two fits if we subtract from the distance modulus of the two models the 
distance modulus /x Mllne (z) of a linearly expanding (but globally flat) homogeneous Milne 
geometry (see e.g [31]), namely if we consider the quantity 

A{z) = /i(z) - fi M[lne {z) = 5 log 

where H is given in units of Mpc -1 . It is clear that positive or negative values of A 
correspond to luminosity distances which are - at a given fixed z - respectively larger 
or smaller than the reference values of the Milne model. 

The case A < is typical of a decelerated Universe like that described by the 
standard cosmological scenario, where, at the same fixed z, the distances are smaller 
(or the received fluxes of radiation, i.e. the apparent magnitudes, are larger) than 
predicted by a linearly expanding model. The case A > 0, on the contrary, corresponds 
at the same z to larger distances (or smaller radiation fluxes) than predicted by linear 
expansion, and is only possible if the model undergoes a period of "effective" accelerated 
expansion. In this last case, the transition across the value A = defines an epoch - 
characterized by the parameter z acc such that A(z acc ) = - marking the beginning of 
the cosmological phase directly imprinted by the kinematic effects of the acceleration. 

The plot of A(z) is presented in Fig. 2 for three cases: the standard CDM- 
dominated (always decelerated) model, and the two best-fit models of Fig. 1 



d L {z) 
IMpc 



5 log 



10 



z(2 + z) 
2H Mpc 



(14) 
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(corresponding to our example of inhomogeneous geometry and to a typical example 
of homogeneous concordance cosmology). In the last two cases we have plotted the 
central values of the fit (solid curves), as well as the corresponding error bands* at the 
95% level of confidence (bounded by the dotted curves). 




Figure 2. The parameter A(z) of Eq. (14) for the two best-fit models of Fig. 1. In 
both cases we have shown the region allowed by the fit at the 95% C.L. (bounded by 
dotted lines). We have also reported (for comparison, and without error band) the 
case of the standard CDM model with O m = 1. 

We can see from Fig. 2 that A(z) is always negative for the CDM model, as 
expected. For the other two models, instead, we have A(z) > in the redshift range 
z < z acc (because, as expected, a successful fit of the SNe data requires the presence of 
a phase describing - or mimicking - accelerated expansion). However, the values of z acc 
defined by the condition A(z acc ) = are largely different in the two models. We find, 
in particular, 

zZ B = 1-07 ± 0.06, 4£D m = L43 ± 10; (15) 

and this difference falls outside the error bands illustrated in Fig. 2 (it is also much 
larger than the experimental uncertainty affecting present redshift measurements). This 
suggests that a precise (near-future?) determination of this parameter could provide 
a clear physical discrimination among different models implementing successful (and 
statistically equivalent) fits of SNe data. 

It should be mentioned, at this point, that in the computations of the error bands we 
have neglected the dispersion of data due to the possible presence of a cosmic background 
of stochastic perturbations: indeed, such a background may induce large errors at very 
small z, but in the range z ~ 1 (typical of z acc ) the induced errors are typically lying 
in the few-percent range [7], hence are not expected to have a crucial impact on the 

* The error bands have been numerically computed by varying the fit parameters within the range 
determined by the corresponding estimated errors. All possible fitting curves lying inside the given error 
region satisfy the constraint: x 2 < x 2 + 5, where x 2 is the value obtained by minimizing Eq.(ll) and 
the constant S depends on the number of parameters and on the confidence level of the fit determined 
by the MINUIT package [29]. In our case, in particular, S = 5.99 for a confidence level of 95%. 
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results illustrated in Fig. 2. The same is expected to be true for the systematic errors 
- possibly slightly bigger than the previous ones, but in any case < 10% - induced on 
z a™ (but not on £acc DM ) by methods of SNe data reduction based on the assumption 
of standard homogeneous cosmology (and used in particular for the Union2 catalogue, 
see e.g. [32]). Finally, we should note that a value of z acc compatible with that of 
the inhomogeneous model considered here could be reproduced also in a homogeneous 
ACDM context, with realistic values of Q m and Q\, but only at the price of introducing 
a large enough negative spatial curvature, with Q k ~ 0.1 (for instance, a model with 
tt m = 0.3, tt A = 0.6, tt k = 0.1 gives ^ DM = 1.087). 

In order to stress the importance of the parameter z acc let us now consider another 
possible form of the phenomenological profile H$(r) appearing in the LTB solution (4), 
for instance the profile* 



H (r) = H + AHtanh 



r — r 
2Ar 



(16) 



We can then explicitly check that different models are characterized by largely different 
values of z acc even within the same class of inhomogenous geometries. By imposing, as 
before, the phenomenological constraint H (0) = 70 Km s -1 Mpc -1 (in order to eliminate 
H), we find that the new profile (16) provides indeed a satisfactory three-parameter fit 
of the Union2 data (see Fig. 3, left panel), with best fit values tq = 2500 ± 322 Mpc, 
Ar = 2387 ± 170 Mpc, AH = 37.5 ± 2.8Kms" 1 Mpc" 1 , at a confidence level of 95%, 
with x 2 /d.o.f. = 1.31. However, the corresponding value of z acc for this model (called 
LTBi in Fig. 3) is significantly different from that of the previous LTB model, and, 
most important, the behaviour of A(z) is exactly the opposite of the standard one, for 
the range of z of our interest (see Fig. 3, right panel). We have checked that the value 
of A(z), for LTB 1; turns back to the standard negative range only for z > 50. 




r =(2500+322) Mpc 
AH=(37.5 ± 2.8) km s ' Mpc" 
Ar=(2387 ± 170)Mpc 



0.5 




1.2 1.4 1.6 



Figure 3. The left panel illustrates the best-fit result for the model characterized by 
the hyperbolic profile of Eq. (16). The corresponding behaviour of A(z), for the range 
of z of interest for this paper, is represented by the curve labelled LTBi reported in 
the right panel. 

* We thank an anonymous referee for this suggestion. See also [33] for other similar profiles. 
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Let us finally comment on the possibility that an off-center position of the 
observer embedded in a spherically symmetric LTB geometry may significantly affect the 
determination of z a ™, thus providing obstructions to a precise discrimination between 
LTB-based and a more conventional (homogeneous) fit of the SNe data. Indeed, if the 
observer is located at a distance r 7^ from the center of a spherically symmetry 
geometry, the corresponding luminosity distance di (referred to the position r ) is no 
longer isotropic but acquires an angular dependence, and this in turn induces an angular 
dispersion of the value of z acc which depends on r , and which obviously grows (in 
modulo) with the growth of r Q . 

The luminosity distance of a source for off-center observers in a LTB geometry has 
been computed in [4] (see also [34]) as a function of z, of the distance tq from the centre, 
and of the polar observation angle 7 (referred to r ). We have applied the results of [4] 
to compute the directional variation of z acc , at fixed values of r®. We have considered, 
in particular, possible displacements from the centre in the range r ^ 10~ 2 r v , because 
- as discussed in [4] - higher values of r would induce a dipole anisotropy too high to 
be compatible with present CMB observations. 

The results of our exercise are illustrated in Fig. 4, where we have plotted the 
fractional variation Az acc /z acc = [-2 acc (r ,7) - z acc (0)]/ z acc (0), for different values of r 
up to 1CT 2 7V, for the LTB model characterized by the parameter z a ™ °f Eq. (15). For 
the normalization of /z Mllne we have consistently used Ho(ro), but we have checked that 
using the fixed value H = 70Kms _1 Mpc _1 simply rescales the zero of the difference 
Az acc , without affecting the overall amplitude of the dispersion. As shown in Fig. 4, the 
angular variation of z acc induced by r 7^ is bounded to be at most at the one-percent 
level, and has thus a negligible impact on the results of Fig. 2. 

0.005 
0.000 

o 

n S -0.005 
< 

-0.010 

* E 3n n 

424 

r 

Figure 4. The fractional variation of the parameter z acc as a function of the angular 
direction 7, for different values of the observer's position r ranging from to 10~ 2 ry. 
The numerical labels of the curves are referred to the values of r , given in units of 
lO- 2 r v . 

In conclusion, we would like to stress again that the inhomogeneous model discussed 
in this paper should not be intended as a realistic alternative to the successful 
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concordance cosmology, but only as a pedagogical example to learn how to distinguish 
different fits of SNe data based on different geometrical schemes. To this purpose we 
have shown, in particular, that in the model of this paper the Universe enters the regime 
directly affected the accelerated kinematics later than predicted by the ACDM scenario, 
i.e. z^™ < -2acc DM - Hence, a precise determination of the transition epoch z &cc (possibly 
through future extensions of the Hubble diagram to higher values of z, or through 
indirect studies of the transfer function of primordial perturbations [35]), could help us 
to physically discriminate among statistically equivalent fits. 
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